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We use the diagrammatic T-matrix approach to analyze the three-body scattering problem between two iden- 
tical fermions and a third particle (which could be a different species of fermion or a boson). We calculate 
the s-wave dimer-atom scattering length for all mass ratios, and our results exactly match the results of Petrov. 
In particular, we list the exact dimer-atom scattering lengths for all available two-species Fermi-Fermi and 
Bose-Fermi mixtures. In addition, unlike that of the equal-mass particles case where the three-body scattering 
T-matrix decays monotonically as a function of the outgoing momentum, we show that, after an initial rapid 
drop, this function changes sign and becomes negative at large momenta and then decays slowly to zero when 
the mass ratio of the fermions to the third particle is higher than a critical value (around 6.5). As the mass ratio 
gets higher, modulations of the T-matrix become more apparent with multiple sign changes, related to the "fall 
of a particle to the center" phenomenon and to the emergence of three-body Efimov bound states. 

PACS numbers: 03.75.Ss, 03.75.Hh, 05.30.Fk 



I. INTRODUCTION 

The dimer-atom scattering process was first solved by Sko- 
rniakov and Ter-Martirosian in 1956 QJJ] , in the context of 
three-nucleon scattering, i.e. a two-body bound state between 
a neutron and a proton (called a deuteron) is scattering with a 
neutron. They considered equal-mass particles, and found that 
while the Born approximation gives a da = 8oab/3 for the 
scattering length between an AB dimer and an A atom, where 
clab is the two-body scattering length between A and B par- 
ticles, the exact result turns out to be clda ~ 118oab JUl- 
This problem has recently been analyzed with the diagram- 
matic T-matrix approach JHHl], in the context of cold quan- 
tum gases (see also J3l in a different context), and the result is 
in perfect agreement with the earlier results. In the same con- 
text, the three-body scattering problem has also been general- 
ized to unequal-mass particles, and analyzed both in real space 
through solving the three-body Schrodinger equation [5J and 
in momentum space through the diagrammatic T-matrix ap- 
proach 1 6], with perfect agreement in between for all mass 
ratios. 

When the mass ratio of the heavy particles to the light one 
is below 13.61, it is also well-established that the three-body 
scattering problem is universal, and that the dimer-atom scat- 
tering length is proportional to the two-body scattering length 
o-ab (just like the equal mass case) with the proportionality 
factor depending only on the masses of the constituent par- 
ticles [5]. However, since three-body Efimov bound states 
emerges for larger mass ratios 17|-|9j], this problem is not uni- 
versal, and an additional parameter coming from the short- 
range (or large-momentum) three-body physics is needed for 
an accurate description. 

Here, we use the diagrammatic T-matrix approach to an- 
alyze the three-body scattering problem between two identi- 
cal fermions and a third particle. Our results are relevant to 
the quantum phases of two-species Fermi-Fermi |[Tol - [l5ll . and 
Bose-Fermi Il6l4l9ll mixtures of atomic gases in the molecu- 
lar limit. The remaining paper is organized as follows. After 
deriving the Skorniakov and Ter-Martirosian integral equation 
generalized for unequal-mass particles in Sec. Ill Al we nu- 



merically solve the resultant equation for all mass ratios in 
Sec. Ill Bl and list the exact dimer-atom scattering lengths for 
all available two-species Fermi-Fermi and Bose-Fermi mix- 
tures in Sec. Ill CI A brief summary of our conclusions is given 
inSecUnl 



II. THREE-BODY PROBLEM 

In this paper, we are interested in the three-body (dimer- 
atom) s-wave scattering between two identical fermions (refer 
to it as A-type particles) and a third particle (refer to it as a 
B-type particle which could be a different species of fermion 
or a boson). In particular, we consider a zero-ranged attractive 
interaction between A and B particles, and assume there is a 
weakly bound resonance between them with the binding en- 
ergy ef, < 0, so that we want to study the scattering between 
this bound state (refer to it as an AB dimer) and the remaining 
A atom. 



A. Dimer-atom scattering T-matrix 

Detailed description of the diagrammatic T-matrix ap- 
proach for the dimer-atom scattering process can be found in 
the literature for equal-mass particles |0-0], and here we give 
details of our calculation [6] for the case of unequal-mass par- 
ticles. 

We begin our analysis by describing the zero temperature 
diagrammatic representation for the dimer-atom scattering T- 
matrix T^ )A (p,po) as illustrated in Fig.Q] where k and p are 
the incoming and outgoing momenta, and oja = k 2 /(2mA) 
and ujd — k 2 /(2mr)) are the frequencies for the incom- 
ing A particle and AB dimer, respectively. Here, we set the 
center-of-mass momentum of the dimer-atom system to zero, 
and e = lja + ^>b + £fc is the total frequency where = 
— lKmABOAs) < is the binding energy of the two-body 
bound state between A and B particles, and mj) = tua + vhb 
and rriAB = 2mBmA/m,D are masses of the AB dimer and 
twice the reduced mass of the A and B particles, respectively. 
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FIG. 1 : Diagrammatic representation of the integral equation for the 
three-body scattering T-matrix Tjf A (p,po) between an AB dimer 
(dimer D consists of one fermionic particle A and a second particle 
B) and a fermionic particle A. Particle B could be a different species 
of fermion or a boson. 

In Fig. Q] single lines represent retarded free propagators 
for the A and B particles 



GU,s(k,w) = — 



(1) 



where loa,b = k 2 / (2mA. b) is the energy, and [ia,b is the 
chemical potential of the corresponding particle. Similarly, 
double lines represent the retarded propagator for the AB 
dimer, and in three dimensions it can be approximated by a 
simple one-pole structure 
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(2) 

which reflects the presence of a two-body bound state between 
A and B particles. 

This dimer propagator is obtained from the resummation of 
the AB polarization bubbles leading to Gc(k, uS) = —g/[l + 
fl , r j 4s(k, ia)], where g > is the strength of the bare interac- 
tion between A and B particles, and the AB polarization bub- 
ble is r A B(k,w) = J2q, qo G A (k + c[,uj + qo)GB{-q,~q ). 
At zero temperature, J2q qo = * f dqdqo / (2ir) 4 in three 
dimensions. Integration over the internal momentum q 
and frequency q leads to Tab^, w) = Tab(O 7 0) + 



3/2 , 
m AB ("D 



1/2 

lo — [i a — (J-b — I (47r), which in com- 

bination with the definition of the two-body scattering length 
&ab = m A bT A B (0,0) /(4tt), and the two-body scattering 
T-matrix T AB (0, 0) = -g/[l + gT A B(0, 0)], between A and 
B particles, give the final result described in Eq. (O. 

In the following, we set ha = Hb = since we are inter- 
ested in the dimer-atom scattering in vacuum. However, note 
that our diagrammatic calculation for the scattering parame- 
ters of the three-body problem is exact, and they are sufficient 
to describe ultracold quantum gases, since experiments are al- 
ways performed at low densities. The calculation of scattering 
parameters in the presence of many other particles (arbitrary 
density) is much more difficult, and it is not discussed here. 
In addition, on the right hand side of Fig. Q] note that the first 
diagram represents a fermion exchange process, i.e. first the 
particle B breaks up with the particle A and then it forms a 
new AB dimer with the remaining A particle. This is the sim- 
plest process contributing to the dimer-atom scattering, e.g. 
Born approximation, and all other (infinitely many) possible 
processes are included in the second diagram. 

In analytical form, the dimer-atom T-matrix T^ A (p,po) 
satisfies the following integral equation 



T£ A (p,p ) = -G B (k + p,u) D -cj A + e b +p ) 

~y^G D (q, Up + e b + qo)G A (~<\,u A - qo) x (3) 
q,<7o 

T D A (q, q )G B (p + q.,u D - u A + tb + Po + qo)- 



The minus signs on the right hand side is due to Fermi-Dirac 
statistics, i.e. exchanging a fermion brings a minus sign unlike 
that of a boson. The integration over frequency qo can be eas- 
ily performed by closing the integration contour in the upper 
half-plane, where both T k D " 4 (q, qo) and Gd(<1,^d + £b + <Zo) 
are analytic functions of qo, and only a simple pole contri- 
bution comes from Ga(— q, u>a ~ qo)- Note that this prop- 
erty of T^ A (q, qo) is due to the form of Eq. Q itself. This 
integration sets qo — (k 2 — q 2 )/ (2mA), and we set po = 
(k 2 — p 2 )j (2mA) in order to have the same frequency depen- 
dence for the T-matrix on both sides of Eq. (01. This leads 
to 
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where we redefine the T-matrix T^^p) = T^^p, (k 2 — p 2 )/(2mA)), and e = k 2 /(2moA) + £fc is the total energy and 
mpA = 2mnmA/ '( m £> + m A) is twice the reduced mass of an AB dimer and an A particle. Since we are interested in the 
zero-ranged low-energy s-wave scattering, we first average out directions of the incoming momentum k, and then of the outgoing 
momentum p, leading to 
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for three-dimensional systems, where T® A (p) = J dfl p J dVL-^T-^ A (p) / (47r) 2 is the angular-averaged T-matrix 
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To obtain the three-body (dimer-atom) s-wave scattering length anA, the total energy e should be set to the binding 
energy et, of the two-body bound state, in the limit of vanishing incoming and outgoing momentum and frequency, i.e. 
auA — 2 , mDAT f ^ A (0)/m AB . This motivates us to introduce the dimer-atom scattering function a^ A {p) for which the in- 
tegral equation becomes 
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in the limit when k — » 0. Here, the dimer-atom scattering 
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gives the full momentum dependence of the dimer-atom scat- 
tering function. 

The integral equation shown in Eq. © as well as the scat- 
tering function expression shown in Eq. © reduce to the 
known results for the equal-mass particles case [2, 3] when 
rriA = ins = m. Since only the fermion exchange pro- 
cess is taken into account in the Born approximation, and that 
neglecting the second term on the right hand side of Eq. (0 
leads to an a = a & ab I m ar, which is consistent with 
the many-body results II20I1 . However, we need to include both 
terms and solve the integral equation in order to find the exact 
dimer-atom scattering length. 



B. Dimer-atom scattering function 
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2: The exact scattering length an a = (0) between an AB 



dimer (dimer D consists of one fermionic particle A and a second 
particle B) and a fermionic particle A is shown as a function of mass 
ratio mA/mB (solid red line) 12 lh . Particle B could be a different 
species of fermion or a boson. Note that the disagreement between 
the exact value and the Born approximation one (dashed line) in- 
creases rapidly as the mass ratio increases. 



Next, we solve numerically the integral equation given in 
Eq. © as a function of the mass ratio mA/mB of the con- 
stituent particles of the dimer 02111 . For this purpose, it is con- 
venient to change the upper integration limit to a finite value 
by a change of variables, e.g. paAB = (1 — x)/(l + %) and 
qa,AB = (1 — y)/(l + y) where 1 > {x,y} > —1. The re- 
sultant integral is calculated by using the Gaussian-Legendre 
quadrature method, and using this discretization, we reduce 
the integral equation to a matrix-eigenvalue problem. 

The exact solutions and the Born approximation values of 
o-da = a o A (0) are shown in Fig. [2] When rriA = tub, we 
find an a ~ l-18aAB, which is in agreement with the results 
previously found for equal-mass particles OH!, HI . The scat- 
tering length a da increases (decreases) from this value with 
increasing (decreasing) mass ratio, and Qda — ► a AB in the 
limit of niA/ms — >■ as expected. These results exactly 
match the few -body results of Petrov |Ht]. It is quite remark- 
able that the diagrammatic T-matrix approach exactly recov- 
ers the few-body results for all mass ratios, since the diagram- 
matic approach is performed in momentum space, while the 
few body approach is performed in real space 12211 . Note that 
the Born approximation values for a da are not in agreement 
with the exact values for any mass ratio, and that the disagree- 
ment increases rapidly with increasing mass ratio, but the gen- 
eral qualitative trend is captured by the Born approximation as 



can be seen in Fig. [2] 

In Fig. [3j we show the scattering function a^ A (p) as a 
function of the outgoing momentum p for some mass ratios 
mA/mB- When tua = tub, this is a monotonically de- 
creasing (positive) function of p with a long tail. However, 
beyond some critical mass ratio, this behavior changes dra- 
matically. For instance, when itla = 10ms, after an initial 
rapid drop, a® A (p) function changes sign and becomes neg- 
ative at paAB ~ 2.72. Beyond this critical momentum, it 
reaches a minimum value and slowly decays to zero from the 
negative side. As the mass ratio gets lower (higher), a^ A {p) 
changes sign at higher (lower) momenta. For instance, when 
rtiA — 6.64tob, corresponding to a three-body system con- 
sisting of one 6 Li atom and two 40 K atoms, it changes sign at 
least once at paAB ~ 14.7 (but it changes sign at least once 
atpaAB ~ 1-96 when tua — 13.61ms). We find clear sign 
changes when tjia > 6.5ms, but we donot know whether 
a o A (p) changes sign for even lower mass ratios since preci- 
sion issues obscure the results. 

Although our calculation is not reliable for mass ratios 
above 13.61 l2lll . for illustration purposes, in Fig. Ob), we 
show ln[a^ ) ' 4 (p)] as a function of p when wia — 30tob. In 
this case, the scattering function first changes sign and be- 
comes negative at paAB ~ 1.25, but then it changes sign again 
and becomes positive at paAB ~ 6.2. Since the large-p behav- 
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FIG. 3: The scattering function a® (p) between an AB dimer (dimer 
D consists of one fermionic particle A and a second particle B) and 
a fermionic particle A is shown as a function of the outgoing mo- 
mentum p for some mass ratios mA/mB 12 ill . Particle B could be 
a different species of fermion or a boson. Here, all length scales are 
in units of the two-body scattering length between A and B particles 
qab, and in (b) we show ln[a,Q A (p)} . 



ior of this function is again obscured due to precision issues, 
as can be seen in the figure, we could not resolve whether it 
has more sign changes at larger p values. However, we men- 
tion that as the mass ratio increases further to 100, aQ A (p) 
changes sign at least 4 times, and the oscillating pattern be- 
comes apparent. 

This effect could be related to the "fall of a particle to the 
center" phenomenon JsJ HHl, and to the emergence of three- 
body Efimov bound states ||7l Hh, which are known to occur 
when the mass ratio is large. When the heavy A particles are 
separated from each other by a distance r <C clab (or equiv- 
alently paAB 3> 1), it is known that an exchange of a light 
B particle mediates an effective —C / (m^r 2 ) attraction be- 
tween A particles. Here, the coefficient C xs 0.162mA/mB 
increases with increasing mass ratio CD. Since A particles 
are fermions, this attraction competes with the Pauli repulsion 
which manifests itself as a centrifugal £(£ + 1) / '(m^r 2 ) po- 
tential, where £ is the angular momentum. Therefore, when 
the masses are comparable (C ~ 0.162), the Pauli repulsion 
is about one order of magnitude stronger than the mediated 
attraction. When the mass ratio ra^/mj reaches a critical 
value, i.e. when C — £(£ + 1), the effective interaction (medi- 



ated attraction plus centrifugal repulsion) between A particles 
vanishes. For the lowest p-wave angular momentum channel 
1=1, this occurs when itia > 12.33tob- 

Beyond this critical mass ratio, there is a second critical 
mass ratio beyond which the effective attraction between A 
particles is strong enough that the A particles stay in an in- 
finitely small region around each other, i.e. the remaining A 
particle falls to the center of the attraction |f23j]. This sec- 
ond critical mass ratio can be approximated by 7 C = 1/4 = 
C — £(£ + 1) 124J. For the lowest p-wave angular momentum 
channel I = 1, this occurs when wia > 13.85tob- Note that, 
since the B particle is already bound to one of the A particles, 
"fall of a particle to the center" is related to emergence of 
three-body Efimov bound states. More accurate calculations 
show that the latter occurs when itia > 13.61ms |18j,|9|]. 

Having presented the diagrammatic T-matrix approach for 
the three-body (s-wave) dimer-atom scattering, for all mass 
ratios mA/m,B, next we discuss the relevance of our results 
to the quantum phases of two-species Fermi-Fermi and Bose- 
Fermi mixtures of atomic gases at ultracold temperatures. 



C. Ultracold Fermi-Fermi and Bose-Fermi mixtures 

It has been shown that l20l l25ll (see i26ll for experimen- 
tal confirmation), in the strong-attraction or molecular limit, 
two-species Fermi-Fermi mixtures with population imbalance 
can be well-described by effective Bose-Fermi models, where 
fermion-fermion pairs behave as molecular bosons (dimers) 
and interact weakly with each other and with the remain- 
ing unpaired (excess) fermions. These simpler models only 
require accurate scattering lengths between two molecular 
bosons (dimer-dimer, i.e. <1dd), an d between a molecular bo- 
son and an unpaired fermion (dimer-atom, i.e. aoA)- Note 
that the exact three- and four-body results in vacuum are suf- 
ficient to describe ultracold atomic mixtures in the molecular 
limit, since experiments are always performed at low den- 
sities. Several fermionic atoms ( 6 Li, 40 K, 87 Sr 12711 . and 
171 Yb [28]) are being currently investigated, and experimental 
methods for studying two-species Fermi-Fermi mixtures are 
being developed in several groups, e.g. 6 Li- 40 K mixture C3- 
[l5tl . Thus, anticipating future experiments involving various 
other mixtures, in Table U we list the exact dimer-atom scat- 
tering lengths for all possible mixtures. 

Similarly, in the strong-attraction or molecular limit, 
two-species Bose-Fermi mixtures (with more fermions than 
bosons, otherwise see B29I0 can be well-described by effec- 
tive Fermi-Fermi models, where boson-fermion pairs can be 
shown to behave as molecular fermions (dimers) and inter- 
act weakly with the remaining unpaired fermions. These sim- 
pler models again only require accurate scattering length be- 
tween a molecular fermion and an unpaired fermion (dimer- 
atom, i.e. clda)- Experimental methods for studying two- 
species Bose-Fermi mixtures are also being developed in sev- 
eral groups Il6l4l9ll . and in Table Ull we list the exact dimer- 
atom scattering lengths for all possible mixtures. 

As we pointed out above, the three-body scattering problem 
is universal for mass ratios below 13.61, and therefore, the 
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TABLE I: The exact scattering length clda = ao >A (0) between a 
bosonic AB dimer (consisting of one fermionic A atom and one 
fermionic B atom) and a fermionic A atom is shown for a list of 
two-species Fermi-Fermi mixtures. Here, clda is in units of the two- 
body scattering length between A and B particles clab , and the mass 
ratios niA/ms are shown inside the parenthesis 12 ill . 



dimer-atom scattering is proportional to the two-body scatter- 
ing length clab with the proportionality factor depending only 
on the masses of the constituent particles. However, since 
three-body Efimov bound states emerges for larger mass ra- 
tios 011, this problem is not universal (an additional param- 
eter coming from the short-range three-body physics is needed 
for an accurate description), and our analysis does not include 
this non-universal effect [5,9]. In particular, the 6 Li- 87 Sr, 6 Li- 
171 Yb and 6 Li- 173 Yb Fermi-Fermi mixtures and 7 Li- 171 Yb, 
and 7 Li- 173 Yb Bose-Fermi mixtures have mass ratios that are 
above the critical ratio for the emergence of three-body Efi- 
mov bound states. 



III. CONCLUSIONS 

In this paper, we used the diagrammatic T-matrix approach 
to analyze the three-body scattering problem between two 
identical fermions and a third particle. The third particle could 
be a different species of fermion or a boson. We calculated the 
exact s-wave dimer-atom scattering length for all mass ratios, 
and our results exactly match the few-body results of Petrov 
who obtained them by solving the three-body Schrodinger 
equation Jit]. It is quite remarkable that the diagrammatic T- 
matrix approach exactly recovers the few-body results for all 
mass ratios, since the diagrammatic approach is performed in 
momentum space, while the few-body approach is performed 
in real space ll22ll . 

We also discussed the relevance of our results to the quan- 
tum phases of two-species Fermi-Fermi and Bose-Fermi mix- 
tures of atomic gases at ultracold temperatures. In particular, 
in the strong-attraction or molecular limit, these mixtures can 
be well-described by simpler effective models, where paired 
atoms behave as dimers and interact weakly with each other 
and with the remaining unpaired atoms. These effective de- 
scriptions require only the scattering lengths between two 
dimers, and between a dimer and an unpaired atom. Antic- 
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(0.99418) 


(1.00583) 


174 Yb 


1.00067 
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(0.49965) 
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176 Yb 


1.00065 
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1.17292 
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(0.03419) 


(0.22714) 


(0.49396) 
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(0.98292) 



TABLE II: The exact scattering length auA = ciq A (0) between a 
fermionic AB dimer (consisting of one fermionic A atom and one 
bosonic B atom) and a fermionic A atom is shown for a list of two- 
species Bose-Fermi mixtures. Here, aoA is in units of the two-body 
scattering length between A and B particles a^s, and the mass ratios 
rriA/mB are shown inside the parenthesis (2lll . 



ipating future experiments, we listed the exact dimer-atom 
scattering lengths for all available two-species Fermi-Fermi 
and Bose-Fermi mixtures. 

In addition, we showed that, unlike that of the equal-mass 
particles case where the three-body scattering T-matrix de- 
cays monotonically as a function of the outgoing momentum, 
after an initial rapid drop, this function changes sign and be- 
comes negative at large momenta and then decays slowly to 
zero when the mass ratio of the fermions to the third particle 
is higher than a critical value (around 6.5). As the mass ratio 
gets higher, modulations of the T-matrix become more ap- 
parent with multiple sign changes. We argued that this effect 
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could be related to the "fall of a particle to the center" phe- (TUBITAK) for financial support, 
nomenon [9|, l23hand to the emergence of three-body Efimov 
bound states HB]. 
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